In-medium nucleon-nucleon potentials in configuration space 
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Based on the thermodynamic Green function approach two-nucleon correlations in nuclear matter at 
finite temperatures are revisited. To this end, we derive phase equivalent effective r-space potentials 
that include the effect of the Pauli blocking at a given temperature and density. These potentials 
enter into a Schrodinger equation that is the r-space representation of the Galitskii-Feynman equa- 
tion for two nucleons. We explore the analytical structure of the equation in the complex fc-plane 
by means of Jost functions. We find that despite the Mott effect the correlation with deuteron 
quantum numbers are manifested as antibound states, i.e., as zeros of the Jost function on the 
negative imaginary axis of the complex momentum space. The analysis presented here is also suited 
for Coulombic systems. 
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I. INTRODUCTION 

Nuclear matter exhibits quite a rich phase structure 
because of strong correlations between nucleons. The for- 
mation of clusters below a certain density (Mott density) 
and the appearance of superconductivity below a critical 
temperature are consequences of these correlations. An 
understanding of nuclear matter at finite temperatures, 
e.g., the knowledge of the equations of state, is largely 
relevant for a description of heavy-ion reactions and as- 
trophysical objects such as the formation and structure 
of neutron stars. 

To describe correlations requires a treatment that goes 
beyond the simple picture of noninteracting quasiparti- 
cles in a background field. A powerful formalism for a 
systematic treatment of correlations is provided by the 
Dyson equations approach, for a review see Ref. El. 
This approach has been used to derive effective in- 
medium equations that can be solved rigorously with 
few-body techniques [PHlq] . For two-body correlations it 
leads to equations known as Galitskii-Feynman or Bcthe- 
Goldstone equations |T4]-fi8[ . These equations include 
the dominant effects of the medium, i.e. the self-energy 
corrections and the Pauli blocking. Since the residual in- 
teraction depends on the Pauli blocking factors the effec- 
tive potential becomes explicitly momentum dependent. 
Therefore the effective equations so far have been solved 
in momentum space only. 

Here we present a solution of the effective in-medium 
two-body equation in co-ordinate space and focus on the 
analytical structure of the residual interaction between 
quasiparticles. Although, in general, solutions should 
not depend on the chosen representation, it is instruc- 
tive to construct an equivalent effective local potential. 
In addition, the techniques presented here are not only 
suited for the two-nucleon problem, but also applicable 
for Coulombic systems where in the simplest case the ef- 
fective interactions arc given in r-space supplemented by 



a Debye screening. 

The interesting region of the phase diagram is that 
close to the Mott transition which indicates a transition 
from a gas of nucleons to a gas of nuclids. Below the Mott 
transition correlations arise as bound states, i.e. a pole 
can be identified for energies E(fi, T, P) < 0. The energy 
depends on the chemical potential /i, the temperature T 
and the cm. momentum P of the cluster. Above the 
Mott transition bound states "disappear" , however, cor- 
relations are still present and can be found by solving the 
respective scattering problem. We explicitely show that 
for the two-nucleon problem, to begin with, a deuteron- 
like in-medium state can be identified as an anti-bound 
state. To arrive at this result we use the technique of 
Jost functions elaborated for local potentials in [|19|-p2[. 
The application to larger clusters will be postponed to a 
future work. 

In Sec. Owe describe the formalism that defines the in- 
medium equations and the resulting medium- modified ef- 
fective nucleon-nucleon potentials for quasiparticles. We 
also briefly recall the Marchenko inversion method to 
derive phase equivalent local potentials and the exact 
method of obtaining potential resonances. Our results 
are presented in Sec. Ill and our conclusions in Sec. IV. 



II. FORMALISM 

A. Effective in-medium Equations 

Within the framework of thermodynamic Green func- 
tions effective in-medium two-body equations known as 
Galitskii-Feynman or Bethe-Goldstone equations (de- 
pending on details) can be derived H]. The residual 
potential that enters into the effective equations, given 
below, includes Pauli blocking factors. The Schrodinger 
type in-medium equation used here is given by 



(H - z)*(12) + ^(1 - h - / 2 )V 2 (12, 1'2')*(1'2') = 0. 

1'2' 

(1) 

where the bare two-nucleon potential is denoted by V%. 
The Matsubara frequency z\ — ir\/(— i(3) + 2fi where 
\x is the chemical potential and [3 = 1/fceT the inverse 
temperature has been analytically continued, i.e z\ — > 
z [HI. The Fermi function is 



/i = f(ki) = 



1 



e /3(e(fei)-A*) _|_ 1 



(2) 



The single-nucleon density p(p, T) can be calculated in 



the standard manner lb 
e that appears in Eq. {\ 



The quasiparticlc self energy 
solves the one particle Green 

function. For an uncorrelated medium the Hartree-Fock 

approximation results in 
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+ A""(k 1 ), (3) 



where 12 denotes proper antisymmetrization of particles 
1 and 2. Hence the effective Hamiltonian of noninter- 
acting identical quasiparticles that enters into eq. (yj) is 
given by 
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For the moment we are interested only in a region of 
rather low density in the vicinity of the Mott transition 
(~ po/10, where po = 0.16 fm~ 3 is normal nuclear matter 
density). We therefore use an effective mass approxima- 
tion for the quasiparticle energy e; i.e. after evaluation 
of Eq. (|j) we fit the effective mass m e ff via 



£ + AHF <« 
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(5) 



The constant shift A^ F can be absorbed in a redefinition 

a hf 



of the chemical potential p, c s = p, — Ag * . Further to 
simplify our analysis we consider the two-body system 
to rest in the medium. It is well known that the influ- 
ence of the medium fades for larger relative momentum 
respective to the medium fL5| . Upon introducing relative 
and cm. coordinates, p and P = 0, for the two-particle 
system the kinetic energy term reads 
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(6) 



where E = z — E cont is the two-body binding/scattering 
energy in the cm. system; the continuum edge for this 
case P — and effective mass approximation is given by 
-Econt = 2A iF jig] . Using a local bare nucleon nucleon 
potential V^r), where r is the conjugate of p, the effective 
potential becomes energy dependent 



V(E,r) = (l-2f(E))V 2 (r), (7) 

where the Fermi function is now explicitly given by 

1 



f(E) 



3 /3(B/2- Mcff ) + l ■ 



(8) 



The effective in-medium Schrodinger-type equation for 
scattering states then reads 



V 



m e g 



+ V{E,r)~E) tp(E,r) = 0. 



(9) 



Note that this equation holds for E > 0, since E — 
p 2 /m e ff for the eigenvalues. 

The strategy is as follows. We study the question 
of what happens to the cluster (for the time being the 
deuteron) for densities above the Mott density and ex- 
plore the location of possible resonances. One may is to 
solve the equations on the unphysical energy sheet (sec- 
ond sheet); for a textbook treatment see, e.g. [E3|. Alter- 
natively, we way utilize the properties of Jost functions 
as elaborated in J19]-g2| . To do so we first derive a phase 
equivalent local potential to Eq. (0) for each chemical 
potential p and temperature T of in terest. This will be 
explained in the next subsection |IIE| . The next step is to 
construct the appropriate Jost function and explore the 
analytical structure by analytic continuation of r to the 



complex plane as explained in the subsection II C . The 



zeros of the Jost functions in the complex "fc-plane" are 
related to bound states, resonances, or anti-bound states 
that can all be identified. 



B. Inverse Scattering Method 

A way to construct the nucleon-nucleus potential in 
configuration space is to use inverse scattering tech- 
niques. The potential in this case is directly obtained 
from the available scattering information and the bound 
states. In the present work we shall employ the inverse 
scattering method at a specific partial wave £ (fixed-^ in- 
version) which has been the subject of several books and 
monographs (see for example, Refs. |24|-[2q|). Therefore, 
we shall outline here only the main features and formulas 
for convenience. 

Most applications at fixed angular momentum em- 
ploy the so-called exactly solvable classes of potentials 
of Bargmann type pTJ . In this approach the S- function, 
defined as a ratio of Jost functions, 



Si(k) 



ft(-k) 
fi(k) 



is parameterized using the rational form 
c ,, \ tt k + ib n -i-r k + a m 



(10) 



(11) 



where Nf, is the number of bound states and M the num- 
ber of terms needed to reproduce the phase shifts from 
to oo. This parametrization ensures the existence of 
bound states at k — ib n and leads to a special class of po- 
tentials known as Bargmann potentials (see Refs. [ p4|j25| 
for more details). They lead to kernels Ki(r,r') of the 
exact Marchenko |2(| and GePfand-Levitan [£8| integral 
equations from which the potentials can be constructed. 
The method has been successfully applied to construct 
nucleon-nucleon potential by von Geramb and collabora- 
tors (see Refs. |2J|pC| and references therein), as well as 
neutron-deuteron potentials 31 32]. 

The Marchenko integral equation, fully equivalent to 
the Schrodinger equation, reads 

/•OC 

K e (r,r') + Ti(r,r') + / K e {r, s)ft(s,r')ds = . (12) 



The kernel Tg(r,r') of this equation is related to the S- 
matrix Si(k), and thus to experiment, via 

T e (r, r') = — hf ] (kr) [1 - S e (k)} h { + ] (kr')dk 
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where hY" (z) is the Riccati-Hankel function, Nf, is the 
number of bound states, and A n £ are the corresponding 
asymptotic normalization constants (in what follows we 
will assume that only one bound state is present in the 
system). Once the Fourier transformation, Eq. (|l3|), is 
obtained, the potential Vt{r) for the partial wave £ is 
calculated from the relation 



Vt(r) 



dK e (r,r) 
dr 



(14) 



With the choice (Till) and for Nf, = 1 the integration in 
(O) can be readily obtained, 



M„ 



F e (r,r r ) = -» £ n e m w+(a e r)wj(a e r') 

m— 1 

+A i wj' (ibir)wf (ifor'), 



(15) 



where TL m are the coefficients from the residues of the 
integrand for the M u poles of (|Tl]) lying in the upper 
complex fc-plane (excluding the one corresponding to the 
bound state ibi). 

Two aspects of the Marchenko scheme should be em- 
phasized. First, in the absence of bound states the 
method provides us with a unique, shape-independent 
(in the sense that the shape is not prechosen) but t- 
dependent interaction that reproduces the phase-shifts 
at all energies. Second, in the presence of bound states 
the potential is not unique as one requires the knowledge 
of the asymptotic normalization constant which can be 
chosen to have any arbitrary value. This will result in 



a set of equivalent local potentials all reproducing the 
scattering data equally well. However, the choice 



Ai = 



-ib) 



ftiib) 



(16) 



where fe(ib) — dfi(k)/dk provide us with a short-range 
unique potential, see Ref. [£4| for a discussion. We men- 
tion here that once the rational fit ( |ll|) is achieved the 
extraction of fe(k) is obtained from the procedure de- 
scribed by Massen et al. |22f | 

In the present work the potentials obtained are energy 
dependent and therefore the bound state for the medium 
modified interaction is not known. However, the fit ( pd| ) 
for the cases where the phase shifts amount to it at the 
origin - irrespective of the number of terms used - always 
provides a constant pole on the positive imaginary fc-axis 
and when ignored (i.e. when it was considered as a pole 
of fi(—k) and thus spurious), highly repulsive potentials 
behaving like ~ 1/r 2 for small r were obtained. These 
are nothin g el se but the shallow supersymmetric partner 
potentials |B2 , B3[ . In contrast, when interpreted as a true 
binding energy with the corresponding asymptotic nor- 
malization constant being calculated via Eq. (16), the fit 
provided us with a short range potential reproducing the 
same scattering data of the energy dependent effective 
potential, Eq. ([?]). 



C. Analytic Properties of the Jost function 

Extraction of the rational parameterization of the Jost 
function does not always provide us with the physical 
poles as some of them may be an artifact of the fitting 
procedure. To extract the actual potential resonances 
one may apply one of the so-called complex energy meth- 
ods. In Refs. |19|-|22|1 a new method has been suggested in 
obtaining the analytical properties of the Jost function 
in the whole complex fc-plane. This method consist of 
replacing the Schrodinger equation by an equivalent sys- 
tem of linear first-order differential equations. By solving 
this system for real momenta k one obtains the Jost so- 
lutions and therefore the Jost function from which the 
phase shifts can be calculated. Solutions for which the 
Jost function has a zero on the imaginary axis provide 
us the bound and anti-bound states (zeros on the positive 
and negative imaginary axis respectively) , while by locat- 
ing a complex zero k r of the Jost function in the fourth 
quadrant of the momentum plane, one obtains the reso- 
nances of the underlying interaction. 

In order to calculate the Jost function for any value of 
momentum on the complex fc-plane we perform a com- 
plex rotation of the coordinate r 



r = x exp(ic/)) x>0 <p e [0, max ] , 



<tt/2 



in the Schrodinger equation and look for a solution of the 
form 



Mk,r) = - [h ( e +) (kr)F^ +) (k,r) + h ( e \kr)F^\k,r) 

(17) 

where h^(kr) are the Riccati-Hankel functions and 
Fg (k,r) are auxiliary functions fulfilling the first or- 
der differential equations jll|] 

A+)/u-\- h< f\ kr )- 



8rF™(k,r) 
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+h'f\kr)F^ ) (k,r) 



^(r)[M +) (fcr)F, (+) (fc,r) 



(18) 



d r Ft\k,r) = _^M^( r ) [h{ + \kr)F^\k,r) 

+h ( f\kr)F^ ] {k,r) 

From the solution of this system one obtains the regu- 
lar solution 4>i(k,r) The use of the system ( |lg| ) has the 
main advantage that it is exact with rigorous boundary 
conditions. 
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h^\kr)ft{k*) + h[ \kr)fi(k)], (19) 
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and that the function F^ coincides asymptotically with 
the Jost function, viz. 



F$-\kr) — ► fe(k) 



(20) 



r — >oo 

More details on the method can be found in Refs. [tL9H22l 

III. RESULTS 

In the present study we employ the central part of 
the Malfliet-Tjon I+III (MT) potential §4), which is 
given in r-space, as the bare potential V2. The three- 
nucleon binding energy for this potential is 8.595 MeV, 
i.e. very close to the experimental value of 8.48 MeV. The 
medium-modified energy dependent potentials V(E J r) 
for five different densities were constructed using Eq. (0). 
The in-medium parameters for these potentials are given 
in Table § 

TABLE I. The potentials used for T = 10 MeV and the 
resulting poles on the imaginary axis of the fe-plane. The 
plus sign for the pole corresponds to bound and the minus to 
anti-bound states. (For the meaning of symbols see text.) 



V(E,r) 


p[fm- 3 ] 


^off [MeV] 


m cff [MeV] 


Pole ib [fm _1 ] 


MT 


- 


- 


939 


+0.2324i 


V003 


0.003 


-22.50 


925.43 


+0.1132i 


V007 


0.007 


-13.29 


907.99 


+0.0055i 


V009 


0.009 


-10.39 


899.69 


-0.0353i 


V017 


0.017 


-2.488 


868.04 


-0.16804 


V034 


0.034 


+7.909 


809.65 


-0.2640i 




FIG. 1. Energy dependent potentials for T — 10 MeV and 
density n — 0.009 fm - (dashed lines) for energies E = 50, 
40, 30, 20, and 10 MeV. Lower energies deviate stronger from 
the isolated MT potential (solid line). 



These are the effective chemical potential /Lt e ff and the 
effective mass m e ff. For simplicity we have used values 
calculated und used earlier M. In Fig. fi] we present the 
potentials V(E,r) for a nuclear density p = 0.009 fm -3 
and for energies E = 10, 20, 30, 40, and 50 MeV and com- 
pare them to the original MT potential. It is seen that 
the potentials become shallower but eventually at higher 
energies the in-medium effects are minimal and the MT 
is practically recovered. The corresponding phase shifts 
for these potentials are shown in Fig. 0. 
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FIG. 2. Si phase shifts for the potential at five differ- 
ent densities (dashed lines). Higher densities deviate stronger 
from the isolated MT phase-shifts (solid line). For densities 
0.009, 0.017, and 0.034 fm -3 the phase-shifts start from zero 
implying that the potential does not sustain a bound state. 
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FIG. 3. The equivalent local potentials V c g(r) obtained 
via Marchenko inversion method for five densities at T = 1 
MeV (dashed lines). Higher densities deviate stronger from 
the isolated MT potential (solid line). 




r[fm] 

FIG. 4. Comparison of the scattering wave functions 
^(E,r) obtained with the _E-independent Marchenko poten- 
tial (dashed lines) with those of V(E,r) (solid lines) for 
n =0.007 fm" 3 at E = 20 and 50 MeV. 



Note that the phase-shifts for the densities 0.009, 
0.017, and 0.034 fm~ 3 start from zero implying that the 
underlying potential does not sustain a bound state. This 
is in agreement with the results of Ref . |l5| , where a sep- 
arable version of the Paris potential has been used as a 
bare nucleon nucleon potential. 

The energy-dependent potentials as given in Eq. (|7]), 
albeit more realistic than the energy independent ones, 
cannot provide us with the bound state or the resonance 
structure of the interaction. To this end we need to con- 
struct an energy independent equivalent local potential 
that can be analytically continued to complex r — ► re ltp 
values, with < tp < 7r/4 and therefore the analytical 
properties of the Jost function in the whole fc-plane can 
be explored. 

To construct energy independent potentials we employ 
the Marchenko inverse scattering method. As discussed 
in Sec. [0] this method requires the phase-shifts at all 
energies. The quality of these potentials depends on the 
cut-off energy E m&K . For sufficiently high £ max (in our 
case ~1500MeV) the off-shell differences of the equiva- 
lent potentials, i.e. the differences of the scattering func- 
tion in the interior region are very small. The resulting 
potentials are presented in Fig. g. Two essential features 
of these interactions should be noted. The attractive well 
becomes deeper but narrower while at the same time a 
repulsion in the interaction region appears as the density 
of the medium increases. This amounts to an decrease 
in the binding energy until, at the Mott density, the two 
nucleons become unbound. Due to the repulsion in the 
interaction region, although quite small, such a poten- 
tial is expected to give rise to resonances. We search for 
resonances in the fourth quadrant of the /c-plane as ex- 
plained in Sec. II C . We find that the bound state, instead 



of being continued off the imaginary axis, moves along 



the negative imaginary fc-axis, i.e. it beconres an anti- 
bound state. Numerical results are given in Table m where 
ib denotes the respective position of the bound (6 > 0) 
or anti-bound (b < 0) state. 

The quality of the potential obtained via inversion is 
demonstrated in Fig. where the scattering wave func- 
tion generated by the energy independent Marchenko po- 
tential V cS (r) for n =0.007 fm~ 3 and at E = 20 and 
50 MeV, is compared with the respective wave function 
for V(E, r) of Eq. (M). The agreement of the respective 
two wave functions is remarkable. Beyond the interac- 
tion region, the two wave functions are, to all practical 
purposes, identical. Similar results were obtained at all 
energies and potentials. 

The potentials obtained via inversion are given in nu- 
merical fornr. Therefore we follow here the same proce- 
dure as used by Massen et al p2|, viz. parameterize the 
potentials using analytical ansatzes, with parameters ob- 
tained via the MERLIN minimization program |R5[. As 
the potentials are intended for use in AGS-type equa- 
tions, we use for this purpose a sum of Yukawa terms 
and the minimization was terminated once an accuracy 
of better than ~ 10~ 5 was achieved at all points. 



IV. CONCLUSIONS 

In the present work we constructed energy-inde- 
pendent but ^-dependent potentials equivalent to the 
medium modified interaction. The purpose of this en- 
deavor is multifold. First, to construct potentials for in- 
terpretation purposes and to find how the medium mod- 
ifications are manifested in r-space. Secondly, to gen- 
erate the necessary numerical methods that can read- 
ily be used to construct energy-independent potentials 



suitable for three- and four-nucleon in-medium calcula- 
tions. Thirdly, to find the analytical properties of the 
corresponding two-body Jost function - the first step to- 
wards the study of resonances in nuclear medium using 
the method described in section IIC. Last but not least, 
we argue that this method although applied in a nuclear 
physics context, can be utilized in other strongly corre- 
lated Fermi systems. 

Two specific features of the constructed interactions 
should be emphasized here; the first one is the move- 
ment of the attractive well towards the interior region 
while at the same time is becomes narrower as the den- 
sities become smaller. The second is the appearance of 
anti-bound states instead of resonances in first quarter 
as the density increases. This reflects a vanishing imag- 
inary part of the respective two-body spectral function. 
Hence, the two-body cluster (even for b < 0) retains the 
character of a quasi-particle, i.e. the correlation does not 
decay. In a more general approach the self-consistency 
requirement (l8) as well as three-body collisions || lead 
to imaginary parts of the two-body spectral function, and 
therefore to a finite width (life time) of the two-particle 
cluster irrespective of b smaller or larger than zero. 

It remains to be seen how these modifications are man- 
ifested in the three- and four-body systems embedded in 
the nuclear medium especially in nuclear scattering, in 
photo- and electro-processes and in the appearance or 
otherwise of resonances in the relevant clusters. 
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